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ABSTRACT 
We present the single-epoch black hole mass estimators based on the C IV A 1549 broad emission line, using 
the updated sample of the reverberation-mapped AGNs and high-quality UV spectra. By performing multi- 
component spectral fitting analysis, we measure the C IV line widths (FWHMciv and line dispersion, ctciv) and 
the continuum luminosity at 1350 A (L1350) to calibrate the C IV-based mass estimators. By comparing with 
the H/3 reverberation-based masses, we provide new mass estimators with the best-fit relationships, i.e., Mbh oc 



r0.50±0.07 
"1350 



cr^iv and Mbh oc L°3fj'=°°''FWHM°iv^°'*'l The new C IV-based mass estimators show significant 
mass-dependent systematic difference compared to the estimators commonly used in the literature. Using the 
published Sloan Digital Sky Survey QSO catalog, we show that the black hole mass of high-redshift QSOs 
decreases on average by ~ 0.25 dex if our recipe is adopted. 
Subject headings: Methods: statistical - Black hole physics - galaxies: nuclei - galaxies: active 



1. INTRODUCTION 

Understanding the growth history of supermassive black 
holes (SMBHs) is one of the fundamental issues in studies of 
galaxy formation and evolution. The intimate connection be- 
tween SMBHs and host galaxies is evidenced through empiri- 
cal correlations between the masses of SMBHs (Mbh) and the 
overall properties of the host galaxy spheroids (e.g., Magor- 
rian et al. 1998; Ferraresse & Merritt 2000; Gebhardt et al. 
2000). The cosmic evolution of these scaling relationships 
has been investigated in the literature, where a tentative evo- 
lution has been reported utilizing observational approaches 
(e.g., Peng et al. 2006; Woo et al. 2006, 2008; Treu et al. 
2007; Merloni et al. 2010; Bennert et al. 2010, 2011; Hiner 
et al. 2012; Canalizo et al. 2012). In order to provide bet- 
ter empirical constraints on the cosmic growth of SMBHs and 
its connection to galaxy evolution, reliable Mbh estimation at 
low and high redshifts is of paramount importance. 

The Mbh can be determined for Type 1 AGN with the rever- 
beration mapping (RM, Peterson 1993) method or the single- 
epoch (SE, Wandel et al. 1999) method under the virial as- 
sumption: Mbh = /^BLR'^y /G, where G is the gravitational 
constant. The size of the broad-line region (BLR), Rblr, can 
be directly measured from RM analysis (e.g., Peterson et al 
2004; Bentz et al. 2009; Denney et al. 2010; Baith et al. 
2011b; Grier et al. 2012) or indirectly estimated from the 
monochromatic AGN luminosity measured from SE spectra 
based on the empirical BLR size-luminosity relation (Kaspi et 
al. 2000, 2005; Bentz et al. 2006, 2009, 2013). The Une-of- 
sight velocity dispersion, AV, of BLR gas can be measured 
either from the broad emission line width in the rms spec- 
trum (e.g., Peterson et al. 2004) obtained from multi-epoch 
RM data or in the SE spectra (e.g.. Park et al. 2012b), while 
the virial factor, /, is the dimensionless scale factor of or- 
der unity that depends on the geometry and kinematics of the 
BLR. Currently, an ensemble average, (/), is determined em- 
pirically under the assumption that local active and inactive 
galaxies have the same Mbh-ct* relationship (e.g., Onken et 
al. 2004; Woo et al. 2010; Graham et al. 2011; Park et al. 
2012a; Woo et al. 2013) and recalibrated to correct for the 

t Author to whom correspondence should be addressed. 



systematic difference of Une widths in between the SE and 
rms spectra (e.g., Collin et al. 2006; Park et al. 2012b). 

The RM method has been applied to a limited sample 
(^ 50) to date, due to the practical difficulty of the extensive 
photometric and spectroscopic monitoring observations and 
the intrinsic difficulty of tracing the weak variability signal 
across very long time-lags for high-z, high-luminosity QSOs. 
In contrast, the SE method can be applied to any AGN if a sin- 
gle spectrum is available, although this method is subject to 
various random and systematic uncertainties (see, e.g., Vester- 
gaard & Peterson 2006, ColHn et al. 2006; McGill et al. 2008; 
Shen et al. 2008; Denney et al. 2009, 2012; Richards et al. 
2011; Park etal. 2012b). 

In the local universe, the SE mass estimators based on the 
H/3 line are well calibrated against the direct H/3 RM results 
(e.g., McLure & Jarvis 2002; Vestergaard 2002; Vestergaard 
& Peterson 2006; Collin et al. 2006; Park et al. 2012b). 
For AGNs at higher redshift (z > 0.7), rest-frame UV lines, 
i.e., Mg II or C IV, are frequently used for Mbh estimation 
since they are visible in the optical wavelength range. Unfor- 
tunately the kinds of accurate calibration applied to H/3-based 
SE BH masses are difficult to achieve for the mass estimators 
based on the Mg II and C IV lines, since the corresponding 
direct RM results are very few (see Peterson et al. 2005; Met- 
zroth et al. 2006; Kaspi et al. 2007). Instead, SE Mbh based 
on these lines can be calibrated indirectly against either the 
most reliable H/3 RM based masses (e.g., Vestergaard & Pe- 
terson 2006; Wang et al. 2009; Rafiee & Hall 2011a) or the 
best calibrated H/3 SE masses (McGill et al. 2008; Shen & Liu 
2012, SL12 hereafter) under the assumption that the inferred 
Mbh is the same whichever line is used for the estimation. 

While several studies demonstrated the consistency be- 
tween Mg II based and H/3 based masses (e.g., McLure & 
Dunlop 2004; Salviander et al. 2007; McGill et al. 2008; 
Shen et al. 2008; Wang et al. 2009; Rafiee & Hall 2011a; 
SL12), the reliability of utilizing the C IV line is still con- 
troversial, since C IV can be severely affected by non-virial 
motions, i.e., outflows and winds, and strong absorption (e.g, 
Leighly & Moore 2004; Shen et al. 2008; Richards et al. 
201 1 ; Denney 2012). Other related concerns for the C IV line 
include the Baldwin effect, the strong blueshift or asymmetry 
of the fine profile, broad absorption features, and the possible 
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presence of a narrow line component (see Denney 2012 for 
discussions and interpretations of the issues). Several studies 
have reported a poor correlation between C IV and H/3 line 
widths and a large scatter between C IV and H/3 based masses 
(e.g., Baskin & Laor 2005; Netzer et al. 2007; Sulentic et 
al. 2007; SL12; Ho et al. 2012; Trakhtenbrot & Netzer 2012). 
On the other hand, other studies have shown a consistency be- 
tween them and/or suggested additional calibrations for bring- 
ing C IV and H/? based masses further into agreement, (e.g., 
Vestergaai-d & Peterson 2006; Kelly & Bechtold 2007; Diet- 
rich et al. 2009; Greene et al. 2010; Assef et al. 2011; Denney 
2012). 

Given the practical importance of the C IV line, which can 
be observed with optical spectrographs over a wide range of 
redshifts (2 < z < 5), in studying high-z AGNs, it is impor- 
tant and useful to calibrate the C IV based Mbh estimators. 
Vestergaard & Peterson 2006 (VP06 hereafter) have previ- 
ously calibrated C IV mass estimators against H/3 RM masses, 
providing widely used Mbh recipes. Since then, however, the 
H/3 RM sample has been expanded and many of RM masses 
have been updated based on various recent RM campaigns 
(e.g., Bentz et al. 2009; Denney et al. 2010; Barth 2011a,b; 
Grier et al. 2012). At the same time, new UV data be- 
came available for the RM sample, substantially increasing 
the quality and quantity of available UV spectra for the RM 
sample. 

In this paper we present the new calibration of the 
C IV based Mbh estimators utilizing the highest quality UV 
spectra and the most updated RM sample. In section 2 we 
describe the sample of H/3 reverberation mapped AGN hav- 
ing available UV spectra. Section 3 describes our detailed 
spectral analysis of the C IV emission line complex to obtain 
the relevant luminosity and line width measurements neces- 
sary for estimating SE Mbh- We provide the updated SE 
C IV Mbh calibration in Section 4 and conclude with a dis- 
cussion and summary in Section 5. We adopt the following 
cosmological parameters to calculate distances in this work: 
Hq = 70 km S-' Mpc"', fl^ = 0.30, and VIa = 0.70. 

2. SAMPLE AND DATA 

For our analysis, we start with the reverberation mapped 
AGN sample, which is considered as a calibration base with 
reliable mass estimates. To date, there are 47 AGNs, for 
which H/3 reverberation based masses are available (Peterson 
et al. 2004; Denney et al. 2006; Bentz et al. 2009; Den- 
ney et al. 2010; Barth et al. 2011a,b; Grier et al. 2012). Of 
those 47 objects, we selected 39 AGNs^ whose archival UV 
spectra are available from International Ultraviolet Explorer 
(lUE) or Hubble Space Telescope (HST^ data archives^. First, 
we collected all available UV spectra covering the C IV spec- 
tral region from the public archives. If there were multiple 
spectra for a given individual object, either multiple epochs 
taken with the same instrument or from multiple instruments, 
we combined the spectra for each instrument by using a stan- 
dard weighted average method to get the better signal-to-noise 
(S/N) ratio. At the same time, we tried to keep contemporane- 
ity as far as possible. Then we selected the best quality spec- 
tra for each object based on visual inspection and by setting a 
limiting S/N ratio of ~ 10 per pixel, which was measured in 
an emission-line free region of the continuum near 1450A or 

^ Note that we also excluded two objects in the list from Peterson et al. 
2004 (i.e., PG 1211+143 and IC 4329A) due to the unreliable RM measure- 
ments as done in VP06, while we included NGC 4593 since a new RM mass 
became available by Denney et al. (2006). 

' http://archive.stsci.edu/iue/ & http://archive.stsci.edu/hst/ 



1700A (see Denney et al. 2013 for the S/N related issues). 
Among these 39 objects, we excluded four AGNs (i.e., NGC 
3227, NGC 4151, PG 1411h-442,PG 1700h-5 18) because they 
are severely contaminated with absorption features. Other 9 
AGNs (i.e., Mrk 79, Mrk 110, Mrk 142, Mrk 1501, NGC 
4253, NGC 4748, NGC 6814, PG 0844H-349, PG 1617H-175) 
were also excluded due to the low quality and unreliability of 
UV spectra. Thus, our sample contains 26 AGNs. Table [T] 
lists the AGNs in the sample and their properties. Note that 
we adopt the updated virial factor, log/ = 0.71 (Park et al. 
2012a; Woo etal. 2013). 

Compared to the previous sample of VP06, one object, 
Mrk 290, is newly included and seven objects (i.e., 3C 120, 
Mrk 335, NGC 3516, NGC 4051, NGC 4593, NGC 5548, 
PG 2130-1-099) have updated reverberation Mbh (Denney et 
al. 2006, Bentz et al. 2009; Denney et al. 2010; Grier et 
al. 2012). One object, PG 0804H-761, that was excluded by 
VP06, is included since it has a new high-quality UV spec- 
trum from the Cosmic Origins Spectrograph (COS) aboard 
HST. Conti-ary to VP06, NGC 4151 is omitted in this work 
due to the strong absorption features near the line center (see 
Section 3). In summary, 13 AGNs have recent high-quality 
UV spectra from HST COS "* compared to VP06. For the re- 
maining objects, UV spectra were obtained from the Space 
Telescope Imaging Spectrograph (STIS) aboard HST for one 
object, from the Faint Object Spectrograph (EOS) aboard HST 
for eight objects, and from Short-Wavelength Prime (SWP) 
camera aboard lUE for four objects as listed in Table |2] We 
corrected the Galactic extinction using the values of E(B-V) 
from the recalibration of Schlafly & Finkbeiner (201 1) listed 
in the NASA/IPAC Extragalactic Database (NED) and the 
reddening curve of Fitzpatrick (1999). 

3. SPECTRAL MEASUREMENTS 

In order to calibrate the C IV Mbh estimator, we mea- 
sured the line width of C IV and the continuum luminos- 
ity at 1350 A, following the multi-component fitting pro- 
cedure developed by Park et al. (2012b) with a modifica- 
tion for the C IV region. We first fitted a single power-law 
continuum using the typical emission-line-free windows in 
both sides of C IV (i.e., - 1340- 1360 A or - 1430- 1470 
A and ~ 1700- 1730 A), which were slightly adjusted for 
each spectrum to avoid the contaminating absorption and 
emission features. We did not subtract the iron emission, 
since it is generally too weak to constrain at least in our 
data sets, although we indeed tested the pseudocontinuum 
model by including the UV Fe II template from Vestergaard 
& Wilkes (2001). After subtracting the best-fit power-law 
continuum, we simultaneously fitted the C IV complex re- 
gion with the multi-component model consisting of a Gaus- 
sian function for the N IV] A1486A, a Gaussian function for 
the Si II A1531Awhenever clearly seen, a Gaussian function 
-I- a sixth-order Gauss-Hermite series for the C IV A 1549 A, 
two Gaussian functions for the He II A1640A, and a Gaus- 
sian function for the O III] A 1663 A. Note that we fitted the 
1600 A feature, which is contaminating the red wing of C IV, 
with a broad He II component (cf. Appendix A. in Fine et 
al. 2010; Marziani et al. 2010). In the fitting process, the 
centers of He II, O III], N IV], and Si II emission line compo- 
nents were fixed to be laboratory wavelengths. We suppressed 
some components in He II, O III], N IV], and Si II lines based 

* For the COS data, we performed a co-addition of multiple exposures with 
the upgraded costools routine (v2.0; Keeney et al. 2012) and then binned the 
spectra to a COS resolution element (~ 0.07 A, ~ 17 km s"') by smoothing 
and re-binning by 7-pixels. 
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Fig. 1 . — Multi-component fitting results for the sample of 26 objects. In eacii panel, tlie observed UV spectrum (black) is overplotted with the best-fit model 
(red), which consists of a single power-law continuum (green), C IV A1549A emission line (magenta). He 11 A1640A emission line (blue), O 111] A1663A emission 
line (brown), and N IV] A1486A emission line (brown). The residual (gray), differences between black and red fines, is presented at tfie bottom of each panel. 



on empirical tests with and without such components. Nar- 
row absorption features were excluded automatically in the 
calculation of x^ statistics by masking out the 3 sigma out- 
liers below the smoothed spectrum (cf. Shen et al. 2011). 
Strong broad absorption features around the C IV line center 
were also masked out manually by setting exclusion windows 
from visual inspection. 

Although it is still controversial whether or not to remove a 
narrow emission-line component from C IVbefore measuring 
the width, we use the full line profile of C IV, i.e., without 
subtracting a narrow emission-line component, in order to be 



consistent with other studies (VP06, Shen et al. 2011, Assef 
et al. 201 1; Ho et al. 2012). We measured the continuum lu- 
minosities at 1350 A and 1450 A from the power-law model 
and measured the C IV line widths (FWHM and a) from the 
best-fit model (i.e., a Gaussian function + a sixth-order Gauss- 
Hermite series) as shown in Figure [T] The measured line 
widths were corrected for the instrumental resolution follow- 
ing the standard practice by subtracting the instrumental reso- 
lution from the measured velocity in quadrature. In Figure 2, 
we explicitly show spectra and best-fit models for the objects 
showing absorption features at the center of C IV. Note that 



Park et al. 



'b 




1500 1550 1600 

Rest Wavelength (A) 

Fig. 2. — The zoom-in view of the C IV for 6 objects, for which fitting 
results are uncertain due to the absorption at the Une center. Black and red 
solid lines represent the observed spectra and the best-fit models, respectively. 



the fitting results are uncertain for these objects, in particular 
NGC4051 (see Section 4.1). 

To assess measurement uncertainties of the line width and 
continuum luminosity, we applied the Monte Carlo flux ran- 
domization method used by Park et al. (2012b; see also Shen 
et al. 2011). Using the 1 , 000 realizations of resampled spec- 
tra made by randomly scattering flux values based on the flux 
errors, we fitted and measured the line width and continuum 
flux, and adopted the standard deviation of the distribution as 
the measurement uncertainties for individual objects as listed 
in Table |2] 

3.1. Continuum Luminosities and Line Widths 
Figure |3] presents the continuum luminosities measured 
at 1350 A and 1450 A, respectively, which are commonly 
adopted for the C IV Mbh estimator Since they are almost 
identical, we choose to use Logo for the mass estimator The 
comparison between FWHM and cr of C IV is plotted in the 
bottom panel of Figure [3] It shows on average, a one-to-one 
relation between FWHM and cr, indicating that the C IV pro- 
file is more peaky than a Gaussian profile, although there is 
large scatter. 

3.2. Comparison to Previous Measurements 

We compare our measurements with those in VP06 in Fig- 
ure IH using the common sample (23 out of 27 objects given 
in their Table 2, except for Mrk 79, Mrk 110, NGC 4151, 
and PG 1617H-175). Since there are multiple measurements in 
VP06, we here show weighted average values of VP06 mea- 
surements for the purpose of comparison. 

For the comparison of Liaso, there is 0.24 dex scatter, which 
may stem from a combination of the differences, e.g., adopted 
spectra and the Galactic extinction correction, between our 
study and VP06. We used the combined single spectra with 
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the best quality while VP06 used all available SE spectra for 
each object. Especially for the objects with the new HST COS 
spectra (red filled circles) observed in different epochs, there 
could be an intrinsic difference due to the variability. In the 
case of the Galactic extinction correction, we utilized the re- 
cent values of E{B-V) listed in the NED taken from the 
Schlafly & Finkbeiner (201 1) recalibration, while VP06 used 
the original values from Schlegel et al. (1998). 

When comparing our ctciv measurements to those of VP06, 
a slight positive systematic trend seems to be present (mid- 
dle panel of Figure HJ. The most likely origin of this trend is 
the difference in the adopted line width measurement meth- 
ods between VP06 and this work. Based on the investigation 
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by Fine et al. (2010), we modeled the C IV complex region 
with multiple components and measured line dispersion from 
the de-blended C IV line model profile, whereas VP06 mea- 
sured line dispersion from the data without functional fits by 
limiting the C IV line profile range to ±10,000 km s~' of the 
line center, regardless of the intrinsic line width of each CIV 
profile. Thus, the line dispersion measured by VP06 will be 
biased if line wings are extended much further than the fixed 
line limit (i.e., underestimation) or the wings are smaller than 
the fixed line limit (i.e., overestimationby including other fea- 
tures). We avoid these biases by de-blending the C IV line 
from other lines using the multi-component fitting analysis 
and measuring the line widths from the best-fit models. Since 
the line dispersion is more sensitive to the line wings than 
the line core, the decomposition and thus recovery of the line 
wing profile from contaminating lines is essential. 

The bottom panel of Figure |4] compares FWHMc iv, indi- 
cating on average consistency between VP06 and this work, 
except for a few outliers. This is because FWHM is less sen- 
sitive to the line wings than ctc iv, hence the difference in the 
measuring method does not generate significant difference in 
measurements. Note that although FWHM is sensitive to the 
nan^ow-line component, both VP06 and our study used the 
full line profile without decomposing the broad and narrow 
components. Instead, another source of discrepancy comes 
from the fact that VP06 measured the line width directly from 
the data while in this study the best-fit models were used for 
line width measurements. Thus, there may be object-specific 
differences depending on how the absorptions above the "half- 
maximum" flux level were dealt with by VP06, and how well 
the functional fits represent the peak of the profile in our study. 

4. UPDATING THE CALIBRATION OF THE C IV SE Mbh 
ESTIMATOR 

By adopting H/3 RM-based masses as true Mbh (see Table 
1), we calibrate the C IV mass estimators by fitting 
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where L1350 is the monochromatic continuum luminosity at 
1350 A and Z\V(C IV) is the line width of C IV, either 
FWHMciv or ctciv. We regress Equation [T] to determine the 
free parameters (q;,/3,7) using the FITEXY estimator imple- 
mented in Park et al. (2012a). Note that this approach is dif- 
ferent from that used by VP06, who adopted the luminosity 
slope from the size-luminosity relation and fixed the velocity 
slope to 2. Instead, this method is consistent with the recent 
approaches described by Wang et al. (2009), Rafiee & Hall 
(2011b), and Shen & Liu (2012). Because a non-linear de- 
pendence is often observed between the line widths of H/3 and 
C IV line (especially based on the FWHM; see Denney 2012 
for a likely physical interpretation), leaving a, /3, and 7 as free 
parameters arguably results in a better statistical regression by 
accommodating the possible covariance between luminosity 
and line width. 

4. 1 . New Calibrations 

In Figure |5] we present the final best-fit calibration results 
for C IV-based mass estimators by directly comparing the 
C IV SE masses with the H/? RM-based masses, using Equa- 
tion[T] The regression results with various conditions and the 
previous calibrations from the literature are listed in Table [5] 
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Fig. 5. — Calibrations of cr-based SE BH masses (top) and FWHM-based 
SE BH masses (bottom) to the H/3 RM-based BH masses. The new sample 
from the recent RM results is marked with a blue open square. The regressed 
parameters (a, 13, j) with the uncertainty estimates are given in the upper left 
part in each panel. 



We adopt the regression results without NGC 4051, which is 
subject to the largest measurement errors among our sample 
since modeling the C IV line of this object is highly uncertain 
due to the strong absorption at the center (see Fig.|2). In addi- 
tion, the variability of NGC 405 1 is expected to have a large 
amplitude since it is the lowest-luminosity object in our sam- 
ple. Thus, NGC 405 1 can add large scatter to the regression 
and potentially skew the slope because there is only a single 
object at the low-mass regime. Thus, excluding this object 
possibly will lead to less biased results in terms of sample se- 
lection and measurement uncertainties. We will present the 
results without NGC 405 1 hereafter unless explicitly stated. 

The slope of the velocity term, when it is treated as a free 
parameter, is closer to the virial assumption (i.e., 2) for the 
crc iv-based estimator, i.e., 7 = 1 .74 ± 0.55 (7 = 1 .45 ± 0.53 
if NGC 405 1 is included) than for the FWHMc iv -based esti- 



mator, i.e., 7 = 0.56 ±0.48 (7 = 0.52 ±0.46 if NGC 4051 is 
included). This reinforces the use of a for characterizing the 
line width of C IV, as suggested by Denney (2012, 2013; see 
also Peterson et al. 2004 and Park et al. 2012b for the case of 
H/3). The slope of the luminosity term (i.e., (3 = 0.51 ± 0.08 
for a; /3 = 0.52 ± 0.09 for FWHM) is almost consistent to that 
of photoionization expectation (i.e., 0.5; Bentz et al. 2006) 
within the uncertainty. This may indicates that asynchronism 
between H/3 and C IV measurements does not introduce a 
significant overall difference for our high-luminosity, high- 
quality calibration sample. 

In this calibration, we treat /3 and 7 as free parameters in ad- 
dition to a. Letting /3 be a free parameter is required to reduce 
luminosity dependent systematics since we are dealing with 
non-contemporaneous Hf3 and C IV measurements, which is 
expected to be not necessarily linear In addition, it is cur- 
rently questionable to directly adopt the C IV size-luminosity 
relation (e.g., Kaspi et al. 2007) for the estimator since it 
is based on such a small sample. Relaxing the constraint of 
7 = 2 for FWHMc iv can be corroborated by the investigation 
by Denney (2012), which shows that there are severe biases in 
measuring FWHMc iv due to the non-variable component and 
dependence on the line shape. These systematic uncertainties 
may be properly calibrated out by taking 7 as a free param- 
eter. In the case of cc iv, however, a similar systematic does 
not seem to be present for ac iv (see Denney 2012). Even if 
we allow 7 to be free, the regression slope for ctc iv is con- 
sistent to the virial expectation (i.e., 2) within la uncertainty, 
thus we opt to fixing 7 to a value of 2, avoiding systematic 
uncertainties due to small number statistics or sample-specific 
systematics. Thus, here we provide the best estimator for the 
C IV-based Mbh (see also, Fig.|5]) as 



log 
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with the statistical scatter against RM masses of 0.29 dex and 
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(3) 



with the statistical scatter against RM masses of 0.35 dex. 
Apart from the interpretation of values of zero point and 
slopes, it is worth noting that these estimators are the best cali- 
brated ones to reproduce Hf3 RM masses as closely as possible 
for the current sample and data sets. 

4.2. Comparison to Previous Recipes 

In Figure |6] we present the systematic difference of Mbh 
estimates based on our new estimators (Equations |2] and |3) 
compared to the previous estimators from VP06 and SL12, 
respectively, using C IV line width and L1350 measurements. 
The ac ly-based Mbh estimates show almost consistent re- 
sults to VP06 with a slight offset of ~ 0.03, which is expected 
from a difference in the adopted values of the virial factor (i.e., 
log/ = 0.71 here versus log/ = 0.74 in VP06). In contrast. 
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the comparison of FWHMc iv-based masses, respectively es- 
timated with our recipe and with that of VP06, shows large 
scatter and a systematic trend that the VP06 recipe underesti- 
mates Mbh in the low-mass regime and overestimates Mbh in 
the high-mass regime, compared to our recipe. The bottom 
panel of Figure |6] shows that the SL12 recipe systematically 
overestimates Mbh over the whole dynamic range of the sam- 
ple (i.e., < 10^ A^o). This is understandable because SL12 
used the FWHMe^g-based Mbh in VP06 as a fiducial one and 
recalibrate the estimator using their high-mass (> 10^ A^o) 
QSOs. Thus, the calibration performed by SL12 in their lim- 
ited dynamical range inherits the overestimation behavior of 
the VP06 recipe with respect to our recipe, and propagates 
it into the low mass regime with larger effect. In addition, 
SL12 subtracted a narrow C IV component before measur- 
ing FWHM of C IV, leading to an overestimated FWHMc iv, 
compared to VP06 and our methods. We note that a large 
dynamic range is necessary for better calibration and investi- 
gation of the biases, as pointed out by SL12. 

In order to explicitly compare the calibration methods used 
in here and VP06, we regress Equation [Hby fixing /3 and/or 7 
with the adopted values in VP06 as listed in Table [3] For the 
(T-based mass estimator, we obtain almost same calibration re- 
sult {a = 6.72 ± 0.06) to that of VP06 (a = 6.73 ± 0.01) using 
the sample including NGC 405 1 . When NGC 405 1 is exclude, 
the zero points reduces slightly (a = 6.69 ± 0.06) and intrinsic 
scatter becomes smaller It is interesting to see the consis- 
tency of the cr -based calibration between our study and VP06, 
despite the syste matic bias in cciv measurements of VP06 as 
shown in Section [3721 We interpret this as follows. Although 
there is a bias in the VP06 measurement method for a, due 
to their choice of line limits (i.e., ±10,000 km s"'), their a- 
based Mbh measurements serendipitously scatter evenly be- 
low and above the central point of the mass scale of the sam- 
ple, consequently resulting in a similar zero point regardless 
of the bias in a measurements. In the end the calibrations are 
very similar, however, the intrinsic scatter of our calibration is 
smaller than that of VP06, which demonstrates a general in- 
crease in accuracy of our cr-based masses over those of VP06, 
advocating for our a measurement prescription. 
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Fig. 6. — Comparison of cr-based (top) and FWHM-based (middle) Mbh 
respectively calculated using the estimators in this study and the estimators of 
VP06. In the bottom panel, Mbh calculated with the FWHM-based estimator 
from SL12 is compared to our mass estimates. Average offset and la scatter 
ai'e given in the lower right comer in each panel. 



4.3. Difference in SMBH Population using the SDSS DR7 
Quasar Catalog 

To demonstrate the effect of our new estimators on Mbh 
studies, we present in Figure |7] the Mbh distribution of the 
SDSS QSO sample as a function of redshift, based on vari- 
ous mass estimators. These masses are calculated using the 
FWHMc IV measurements from Shen et al. (2011), who pro- 
vides only FWHM measurements using SDSS DR7 spectra. 
Note that FWHMc iv -based Mbh determined with our new 
estimator is on average smaller by ^ 0.25 dex than that cal- 
culated with the previous estimator by VP06, since the VP06 
recipe tends to overestimate Mbh in the high-mass regime as 
explained in Section 4.2. In contrast, there is a smooth tran- 
sition between Mg Il-based masses estimated from the recipe 
of Wang et al. (2009) and C IV-based masses from our new 
calibration since both estimators are based on the same cali- 
bration scheme (see Section|4]i. 

Kelly & Shen (2013) derived a predicted maximum Mbh 
of broad-line QSOs as a function of redshift (their Figure 7), 
showing a slight trend that the maximum Mbh was larger at 
higher redshift. However, this subtle trend may simply be 
a result of systematic overestimation of C IV masses at high 
redshift because their mass determination was based on the 
VP06 recipe. Adopting our new C IV mass estimator may 
eliminate such a trend. 
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5. DISCUSSION AND CONCLUSIONS 

We investigated the calibration of C IV Mbh estimators 
based on the updated sample of 26 AGNs, for which both 
H/3 reverberation masses and UV archival spectra were avail- 
able. The sample of AGNs with RM masses as well as the 
UV spectra including C IV have been expanded and updated 
since the calibrations performed by VP06; it is therefore use- 
ful to revisit the calibration of the C IV-based mass estima- 
tors to provide the most consistent virial Mbh estimates using 
the C IV emission line. Major differences of the calibration 
method between VP06 and the current study is twofold. First, 
we derived line widths (i.e., line dispersion and FWHM) from 
the spectral fits by performing multi-component fitting on the 
C IV complex region to accurately de-blend C IV from other 
contaminating lines while VP06 measured the line width of 
CIV directly from the spectra. When "applying" a SE scaling 
relationship to calculate Mbh, it is important to use the same 
fitting and line width measurement prescriptions that were 
used in "calibrating" the scaling relationship because signifi- 
cant systematic differences can arise in Mbh estimates if dif- 
ferent analysis and measurement techniques are utilized (e.g., 
Assefetal. 2011, Denney 2012, SL 12, Park et al. 2012b). 
Second, we treated the slope parameters (i.e., /? and 7) in the 
virial Mbh equation (i.e., Eq. [Til as free parameters as in Wang 
et al. (2009), which is particularly important for FWHM. 

We provided the best-fit calibrations for both ctc iv- and 
FWHMc IV-based mass estimators. While we presented a 
consistent estimator for the cr-based masses to that of VP06, 
we obtained significantly different Mbh estimator for the 
FWHM-based masses, presumably due to relaxing the con- 
straint of the virial expectation (i.e., 7 = 2) to mitigate the 
FWHM-dependent biases. We generally recommend to use 
the cr-based mass estimator if the a measurement is avail- 



able, as it shows the better consistency with the virial relation, 
and (7 -based masses show a smaller scatter than the FWHM- 
based masses when compared to H/3 RM-based masses. Us- 
ing (Tc IV-based estimator is also preferred by Denney (2012), 
who showed that the ctc iv measured from mean spectra is 
the better tracer of the broad-line velocity field than the 
FWHMc IV since FWHM of C IV is much more affected by 
the non-variable C IV core component. 

Compared to the previously calibrated FWHMc iv estima- 
tor by VP06, our new estimator shows a systematic trend as 
a function of mass. The VP06 recipe overestimates Mbh in 
the high-mass regime (i.e., > 1O**M0) while it underestimate 
Mbh in the low-mass regime (i.e., < IO^Mq), compared to 
C IV masses based on our new estimator This systematic dis- 
crepancy is due to a combination of effects, including differ- 
ence in the RM sample and updated RM masses, newly avail- 
able UV spectra, emission-line fitting method, and calibration 
method. For the SDSS quasar sample (Shen et al. 2011), 
we find that Mbh estimates based on our new estimator are 
systematically smaller by ^ 0.25 dex than those based on the 
previous recipe of VP06. 

One of the main differences in calibrating the FWHMc iv- 
based mass estimator is that we fit the exponent of velocity 
(/3 = 0.56) as in Eq. 3, instead of adopting /3 = 2 as in VP06. 
This provides effectively the same effect as adopting a vary- 
ing virial factor If a constant virial factor is used for mass de- 
termination, then FWHM-based masses will show systematic 
difference compared to cr-based masses, since the FWHM/cr 
ratio has a broad distribution, while the fiducial RM masses 
are derived with a measurements fromrms spectra. To resolve 
this issue, Collin et al. (2006), for example, introduced vary- 
ing virial factors for the FWHMh^ -based mass estimator de- 
pending on the range of the line widths. In our case, relaxing 
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the virial (FWHM ) requirement in calibrating FWHMc iv- 
based masses against o'h/3 -based RM masses provides virtu- 
ally the same effect as adopting a varying virial factor, thus 
resulting in better consistency with cr-based masses. It also 
mitigates the bias caused by the contamination of the non- 
variable C IV emission component, where FWHMc iv-based 
masses in objects with 'peaky' ('boxy') profiles are under- 
cover-) estimated with previous FWHMc ly-based mass esti- 
mators. 

The calibration of mass estimators provided in this study 
still suffers from a sample bias as in the case of VP06. The 
incompleteness or lack of low-mass objects (i.e., < IO^Mq) in 
the current sample will be resolved when new HST STIS ob- 
servations become available for six low-mass reverberation- 
mapped AGNs (GO-12922, PI: Woo). However, the extrap- 
olation of this calibration to high-luminosity, higher-redshift 
AGNs more similar to the SDSS sample can only be realized 
with the extension of the RM sample to this regime — an en- 
deavor that we strongly advocate. 

Apart from the calibration analysis performed in this and 
previous studies, several schemes to correct for the C IV- 
based masses have been suggested in the literature to re- 
duce the large scatter between the C IV-based masses and 
the H/3-based masses. For example, Assef et al. (2011) 
suggested a prescription to reduce mass residuals using the 
ratio of the rest-frame UV to optical continuum luminosi- 
ties based on a sample of 12 lensed quasars. Shen & Liu 
(2012) reported a poor correlation with large scatter be- 
tween FWHMh;3 and FWHMciv using a sample of 60 high- 
luminosity QSOs, showing that some part of the scatter cor- 
related with the blueshift of C IV with respect to H/3. They 
suggested a correction for the FWHM of C IV and C III] lines 
as a function of the C IV blueshift. Recently, Denney (2012) 
showed that the C IV line profile consists of both non-variable 
and variable components based on the sample of seven AGNs 
with C IV reverberation data, and concluded that this non- 
variable component is a main source of the large scatter of 
the C IV SB Mbh- They provided an empirical correction for 



the FWHM-based mass depending on the line shape as param- 
eterized as the ratio of FWHM to the line dispersion. Since 
the C IV line region is more likely to be affected by non-virial 
motions such as outflows and winds than the lower ionization 
line region, such as H/3 (e.g., Shen et al. 2008; Richards et 
al. 2011), aforementioned corrections are also important and 
worth investigating further with a larger sample with enlarged 
dynamic range. 

In general, correcting for possible systematic biases and 
providing accurate Mbh estimates is crucial for studies of 
the cosmic evolution of BH population, particularly at high- 
redshifts (e.g.. Fine et al. 2006; Shen et al. 2008, 2011, 
Shen & Kelly 2012; Kelly & Shen 2013). Thus, it is im- 
portant to ensure a reliable calibration at the high-mass end 
(> lO'' M0) since the C IV mass estimators are most applica- 
ble to high-mass AGNs at high-redshift utilizing optical spec- 
troscopic data from large AGN surveys. Note that the current 
RM sample used for calibrating C IV mass estimators still suf- 
fers from the lack of high-luminosity AGNs, suggesting that 
the RM sample may not best represent the high-luminosity 
QSOs at high-redshifts, i.e., SDSS QSOs. Thus, obtaining di- 
rect C IV reverberation mapping results for high-mass QSOs 
will be even more useful (see Kaspi et al. 2007 for tenta- 
tive results), despite the practical observational challenges. 
Such measurements will be used to better determine a reli- 
able C IV size-luminosity relation and to directly investigate 
non-varying component of C IV. 
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TABLE 1 
Optical spectral properties from H/3 reverberation mapping 
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2539 ±466 


7.80 ±0.04 


6 


3C390.3 
Arkl20 


0.05610 
0.03230 


23.60!|g 


3105 ±81 
1896 ±44 


9958 ±1046 
5364 ±168 


8.43 ±0.10 
8. 14 ±0.06 


1 
1 


Fairall9 


0.04702 


17.40+3^^1 


3787 ±197 


6901 ±707 


8.38 ±0.10 


1 


Mrk279 


0.03045 


16J0^1Z 


1420 ± 96 


3385 ±349 


7.51 ±0.12 


1 


Mrk290 


0.02958 


8.72!!;^^ 


1609 ±47 


4270 ±157 


7.36 ±0.07 


4 


Mrk335 


0.02578 


14 1+0.4 
79.60+|-| 

19.05+1*^ 
11.68! 1 

IQ 20+2.80 

73+8:^^ 


1293 ±64 


1025 ±35 


7.37 ±0.05 


6 


Mrk509 


0.03440 


1276 ±28 


2715 ±101 


8. 12 ±0.04 


1 


Mrk590 


0.02638 


1653 ±40 


2512±101 


7.65 ±0.07 


1 


Mrk817 


0.03145 


1636 ±57 


3992 ± 302 


7.66 ±0.07 


1 


NGC3516 

NGC3783 


0.00884 
0.00973 


1591 ±10 

1753 ±141 


5175 ±96 
3093 ±529 


7.47 ±0.05 
7.44 ±0.08 


4 
1 


NGC4051 
NGC4593 


0.00234 
0.00900 


927 ± 64 
1561 ±55 


1034 ±41 
4141 ±416 


6.20 ±0.14 
6.96 ±0.09 


4 
2 


NGC5548 
NGC7469 


0.01717 
0.01632 


4 lg+0.86 

4 50+^^^^ 
4.3U „„ 

111.00+f2» 


3900 ±266 
1456 ± 207 


12539 ±1927 
2169 ±459 


7.80 ±0.14 
7.05 ± 0.05 


3,5 
1 


PG0026+129 


0.14200 


1773 ±285 


1719 ±495 


8.56 ±0.11 


1 


PG0052+251 


0.15500 


89.80+2«o 
146 90+18.85 

150.10+11 


1783 ±86 


4165 ±381 


8.54 ±0.09 


1 


PG0804+761 
PG0953+414 


0.10000 
0.23410 


1971 ±105 
1306 ±144 


2012 ±845 
3002 ± 398 


8.81 ±0.05 
8.41 ±0.09 


1 

1 


PG 1226+023 


0.15830 


306.80!™:™ 


1777 ±150 


2598 ±299 


8.92 ±0.09 


1 


PG1229+204 
PG 1307+085 
PG1426+015 


0.06301 
0.15500 
0.08647 


37.80+21*! 
105.60+11 
95.00+11; 


1385 ±111 
1820 ±122 
3442 ± 308 


3415 ±320 
5058 ±524 
6323 ±1295 


7.83 ±0.23 
8.61 ±0.12 
9.08 ±0.13 


1 

1 
1 


PG1613+658 


0.12900 


40.io!|5:|o 


2547 ± 342 


7897 ±1792 


8.42 ±0.22 


1 


PG2 130+099 


0.06298 


12.8!i:^ 


1825 ±65 


2097 ± 102 


7.63 ±0.04 


6 



Note. — Col. (1) Name. Col. (2) Redshifts are from the NASA/IPAC Extragalactic Database (NED). 
Col. (3) Rest-frame H/3 time lag measurements. Col. (4) Line dispersion (crYms) measured from rms 
spectra. Col. (5) FWHM measured from rms spectra Col. (6) Mbh estimates from reverberation mapping: 
Mbh(RM) = /cTcentO"jn,s/G where the virial factor / is adopted from Park et al. (2012a) (i.e., log/ = 0.71). 
Col. (7) References. 1. Peterson et al. 2004; 2. Denney et al. 2006; 3. Bentz et al. 2009; 4. Denney et al. 
2010; 5. Park et al. 2012b; 6. Grier et al. 2012. 

" Note that as in VP06 the Mbh taken from Table 8 in Peterson et al. (2004) are based on the weighted 
average virial products from RM results of all reliable different emission lines rather than H/3 only results. 
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TABLE 2 
Ultraviolet spectral properties from C iv single-epoch estimates 



Object 


Telescope/Instrament 


Date Observed 


S/N 


E(B-V) 


log(AL;^/ergs-l) 


o-SE 


FWHMsE 


log(MBH/M0) 


log(MBH/M0) 


Notes 










(1450A or 1700A) 




(1350A) 


(CIV) 


(CIV) 


((T(C IV), SE) 


(FWHM(C IV), SE) 












(pix-') 


(mag) 




(kms-') 


(kms-') 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 




3C120 


lUE/SWP 


1994-02-19,27;1994-03-ll 


12 


0.263 


44.399 ±0.021 


3106 ±157 


3093 ±291 


7.89 ±0.05 


7.96 ± 0.03 






3C390.3 


HST/FOS 


1996-03-31 


18 


0.063 


43.869 ±0.003 


6154 ±65 


5645 ± 202 


8.22 ±0.01 


7.83 ±0.01 






Arkl20 


HST/FOS 


1995-07-29 


17 


0.114 


44.400 ±0.005 


3219±53 


3471 ±108 


7.93 ±0.01 


7.99 ±0.01 






Fairall9 


HST/FOS 


1993-01-22 


24 


0.023 


44.442 ± 0.004 


2694 ± 20 


2649 ± 77 


7.79 ±0.01 


7.95 ±0.01 






Mrk279 


HST/COS 


2011-06-27 


9 


0.014 


43.082 ±0.004 


2973 ± 53 


4093 ± 388 


7.20 ±0.02 


7.35 ±0.02 






Mrk290 


HST/COS 


2009-10-28 


24 


0.014 


43.611 ±0.002 


3531 ±32 


2052 ±36 


7.61 ±0.01 


7.45 ±0.01 






Mrk335 


HST/COS 


2009-10-3 1;2010-02-08 


29 


0.032 


43.953 ±0.001 


1876 ±12 


1772 ±14 


7.23 ±0.01 


7.59 ±0.01 






Mrk509 


HST/COS 


2009-12-10,11 


107 


0.051 


44.675 ±0.001 


3568 ±9 


3872 ±18 


8.15 ±0.01 


8.16±0.01 






Mrk590 


lUE/SWP 


1991-01-14 


17 


0.033 


44.094 ±0.007 


3479 ±165 


5362 ±266 


7.84 ±0.04 


7.94 ±0.01 






Mrk817 


HST/COS 


2009-08-04;2009-12-28 


38 


0.006 


44.326 ±0.001 


3692 ±23 


4580 ±48 


8.01 ±0.01 


8.02 ±0.01 






NGC3516 


HST/COS 


2010-10-04;201 1-01-22 


20 


0.038 


42.615 ±0.002 


4006 ±49 


2658 ± 34 


7.22 ±0.01 


7.00 ±0.01 


abs 




NGC3783 


HST/COS 


2011-05-26 


29 


0.105 


43.400 ±0.001 


2774 ±91 


2656 ±444 


7.30 ±0.03 


7.41 ±0.04 




n 


NGC4051 


HST/COS 


2009-12-11 


23 


0.011 


41. 187 ±0.001 


1290 ±139 


1122 ±309 


5.53 ±0.09 


6.05 ± 0.07 


abs 


p_ 


NGC4593 


HST/STIS 


2002-06-23,24 


10 


0.022 


43.761 ±0.005 


2946 ±162 


2952 ±166 


7.53 ±0.05 


7.62 ±0.01 




S^ 


NGC5548 


HST/COS 


2011-06-16,17 


36 


0.018 


43.822 ±0.001 


4772 ±80 


1785 ±82 


7.98 ±0.01 


7.53 ±0.01 


abs 


p 


NGC7469 


HST/COS 


2010-10-16 


32 


0.061 


43.909 ±0.001 


2849 ± 237 


2725 ± 66 


7.57 ± 0.07 


7.68 ±0.01 




3' 


PG0026+129 


HST/FOS 


1994-11-27 


25 


0.063 


45.236 ±0.005 


4965 ±113 


1604 ±50 


8.72 ±0.02 


8.24 ±0.01 




era 


PG0052-I-251 


HST/FOS 


1993-07-22 


21 


0.042 


45.292 ±0.004 


4648 ± 50 


5380 ±87 


8.69 ±0.01 


8.56 ±0.01 




n 


PG0804-I-761 


HST/COS 


2010-06-12 


34 


0.031 


45.493 ±0.001 


2585 ± 20 


3429 ± 23 


8.28 ±0.01 


8.56 ±0.01 




< 


PG0953-I-414 


HST/FOS 


1991-06-18 


18 


0.012 


45.629 ±0.005 


3448 ±55 


3021 ±74 


8.60 ±0.01 


8.60 ±0.01 




1 

a' 


PG 1226+023 


HST/FOS 


1991-01-14,15 


93 


0.018 


46.309 ±0.001 


3513±29 


3609 ± 29 


8.96 ±0.01 


8.99 ±0.01 




C/5 


PG 1229+204 


lUE/SWF 


1982-05-01,02 


28 


0.024 


44.609 ± 0.009 


2621 ±90 


4023 ±163 


7.85 ±0.03 


8.14±0.01 




p. 


PG 1307+085 


HST/FOS 


1993-07-21 


14 


0.030 


45. 113 ±0.006 


4237 ± 80 


3604±111 


8.52 ±0.02 


8.37 ±0.01 




D3 


PG1426+015 


IUE/SW¥ 


1985-03-01,02 


45 


0.028 


45.263 ±0.004 


4808 ± 305 


4220 ±25 8 


8.71 ±0.06 


8.49 ±0.02 




PG1613+658 


HST/COS 


2010-04-08,09,10 


37 


0.023 


45.488 ±0.001 


4204 ±17 


6398 ±51 


8.70 ±0.01 


8.71 ±0.01 




a 


PG2 130+099 


HST/COS 


2010-10-28 


22 


0.039 


44.447 ±0.001 


2225 ±47 


2147 ±18 


7.63 ±0.02 


7.90 ±0.01 







Note. — Col. (1) Name. Col. (2) Telescope/Instrument from which archival UV spectra were obtained. Note that the new COS spectra were obtained after 2009. Col. (3) Observation date for 
combined spectra. Col. (4) Signal-to-noise ratio per pixel at 1450A or 1700A in rest-frame. Col. (5) E(B-V) are from the NASA/IPAC Extragalactic Database (NED) based on the recalibration of 
Schlafly & Finkbeiner (2011). Col. (6) Continuum luminosity measured at 1350A. Col. (7) Line dispersion ((Time) measured from SE spectra. Col. (8) FWHM measured from SE spectra. Col. (9) SE 
mass estimates based on an^^ from the new estimator. Col. (10) SE mass estimates based on FWHM from the new estimator. Col. (11) abs: The C IV line profile is affected by a strong absorption near the 
line center. Thus the emission-line modeling could be uncertain although that region was masked out. 
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TABLE 3 
C IV Mass Calibration results 
log[MBH(RM)/MQ] = a + (5 logCLijso/lO"'' erg s"') + 7 log[ziV(C IV)/1000 km s" 



AV(C IV) 



mean offset 
(dex) 



la scatter 
(dex) 



Ref. 



Previous calibrations 



""line 

0"line 

FWHM 

FWHM 

FWHM 



6.73 + 0.01 
6.73 + 0.02 
6.66 + 0.01 
6.69 + 0.01 
8.021 



0.53 
0.53 
0.53 
0.53 
0.471 



2 
2 
2 
2 
0.242 



0.33 
0.37 
0.36 
0.43 



This work 



0.03 



0.28 



vpoe" 

VP06'' 

vpoe" 

VP06'' 

SL12 



""line 

FWHM 



7.02 + 0.29 
7.51+0.21 



0.46 + 0.07 
0.51+0.08 



1.45 + 0.53 
0.52 + 0.46 



0.30 + 0.04 
0.34 + 0.05 



0.00 
0.00 



0.298 
0.340 



This work (fixing 7 = 2) 










o-l,„e 6.75 + 0.07 0.43 + 0.07 
FWHM 6.81+0.09 0.40 + 0.08 


2 
2 


0.30 + 0.05 
0.41+0.05 


0.00 
0.00 


0.308 
0.421 


This work (fixing fi = 0.5) 










(Tline 7.07 + 0.27 0.5 
FWHM 7.50 + 0.18 0.5 


1.34 + 0.48 
0.54 + 0.36 


0.29 + 0.04 
0.33 + 0.04 


0.00 
0.00 


0.301 
0.340 


This work (fixing /3 = 0.5 and 7 = 2) 










(Time 6.73 + 0.06 0.5 
FWHM 6.78 + 0.08 0.5 


2 
2 


0.30 + 0.04 
0.41+0.05 


0.00 
0.00 


0.318 
0.439 


This work (fixing fi = 0.53 and 7 = 2 for the comparison 


to VP06) 






(Ti,„e 6.72 + 0.06 0.53 
FWHM 6.77 + 0.08 0.53 


2 
2 


0.31+0.04 
0.42 + 0.05 


0.00 
0.00 


0.328 
0.449 


This work (w/o NGC 4051) 










(Time 6.84 + 0.30 0.51+0.08 
FWHM 7.48 + 0.24 0.52 + 0.09 


1.74 + 0.55 
0.56 + 0.48 


0.29 + 0.04 
0.35 + 0.05 


0.00 
0.00 


0.293 

0.347 best-fit"^ 


This work (w/o NGC 405 1 ; fixing 7 = 2) 










(Time 6.71 + 0.07 0.50 + 0.07 

FWHM 6.78 + 0.10 0.45 + 0.10 


2 

2 


0.28 + 0.04 
0.41+0.05 


0.00 
0.00 


0.295 best-fit"^ 
0.419 


This work (w/o NGC 405 1 ; fixing 13 = 0.5) 








(Time 6.84 + 0.29 0.5 
FWHM 7.47 + 0.23 0.5 


1.75 + 0.51 
0.59 + 0.44 


0.28 + 0.04 
0.34 + 0.04 


0.00 
0.00 


0.293 
0.347 


This work (w/o NGC 405 1 ; fixing fJ = 0.5 and 7 = 2) 








(Time 6.71+0.06 0.5 
FWHM 6.75 + 0.08 0.5 


2 
2 


0.28 + 0.04 
0.40 + 0.05 


0.00 
0.00 


0.295 
0.422 


This work (w/o NGC 405 1 ; fixing /3 = 0.53 and 7 = 2 for the comparison 


to VP06) 




(Time 6.69 + 0.06 0.53 
FWHM 6.74 + 0.08 0.53 


2 
2 


0.28 + 0.04 
0.41+0.05 


0.00 
0.00 


0.296 
0.425 



Note. — The mean offset and Icr scatter for our calibrations are measured from the average and standard 
deviation of mass residuals between RM masses and calibrated SE masses, A = logMBH(RM) — logMBH(SE). 
" The VP06 calibration based on all individual measurements for each object. 
'' The VP06 calibration based on weighted averages of all measurements for each object. 
'^ We suggest these calibrations as the best Mbh estimators. The viral factor is assumed as log/ = 0.71. 



